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Abstract 

We propose an approach to statistical systems on lattices with 
sphere-like topology. Focusing on the Ising model, we consider the 
thermodynamic limit along a sequence of lattices which preserve the 
fixed large scale geometry. The hypothesis of scaling appears to hold 
at criticality, pointing at a sensible definition of the continuum limit 
of the model in the curved space. 
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The elucidation of the effects of gravity on statistical systems has already 
deserved a lot of attention. The discussion has been mainly centered on 
the influence of quantizing the spatial geometry of two-dimensional statisti- 
cal models at criticality. The determination of the gravitational dressing of 
conformal weights can be considered one of the most important theoretical 
breakthroughs of the past decade. Let us recall that it was first achieved 
in the light-cone gauge formulation of quantum gravity [l|], and later by 
use of conformal field theory methods ||. The physical problem which is 
adressed by this study is that of considering cooperative phenomena on a 
fluctuating surface, and it is in this kind of setting that the gravitational 
scaling dimensions have been confirmed by the method of dynamical trian- 
gulation of random surfaces || |, A different topic arises, however, by 
considering the effect on the critical behaviour of a fixed spatial geometry. 
In this respect, it is important to stress that the effects of curvature have 
not been assessed on statistical systems. Since these are made of a discrete 
number of lattice points, the main drawback in posing the problem concerns 
the way in which the thermodynamic limit can be approached keeping fixed 
the spatial geometry. Focusing in two- dimensions, it is not obvious how to 
enforce this constraint, specially when the mean curvature over the surface is 
positive. This corresponds to having a positive Euler number, which forces 
the topology to be that of the sphere. On the other hand, this is the first 
case of interest since quite significant physical systems are already known 
which realize the above condition on the curvature, from a simplicial point of 
view. The most notorious instance is the buckminsterfullerene, the molecule 
with sixty carbon atoms placed at the vertices of a truncated icosahedron 
|||. Recently, carbon aggregates have been discovered in which each shell 
has an increasing number of atoms disposed on a closed lattice of sphere-like 
topology [0]. These kind of systems invite to open the study of statistical 
properties in the presence of nonvanishing curvature. 

From the theoretical point of view, the basic question is to give a growing 
pattern by which a lattice can be built systematically from a previous one, 
keeping the same distribution of the curvature on the new scale. In this letter 
we propose a definition of the thermodynamic limit along some simple classes 
of two-dimenional lattices with the topology of the sphere. According to 
simplicial geometry, one assigns to every lattice face with number of vertices 
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Hi the surface element 

ASi = -rii 
and a value of the curvature given by 



(1) 



The search for the mentioned growing pattern leads to consider hexagonal 
lattices with a finite number of defects. Let us constrain the problem by 
allowing only convex (B4 > 0) dislocations. It turns out that these have to 
be placed according to the spatial orientation of the vertices of the tetrahe- 
dron, the octahedron or the icosahedron. To be more precise, the lattices 
we consider have as building blocks triangular pieces of honeycomb lattice 
of the type shown in Figure 1, which are to be assembled as the faces of the 
given regular polyhedron M. When this is done, we end up with a closed 
lattice with constant coordination number in which all the curvature is con- 
centrated at a finite number of locations. These correspond to four three-fold 
rings in the case of the tetrahedron, six four-fold rings in the octahedron and 
twelve five-fold rings in the icosahedron. Given the polyhedron, though, the 
family of lattices of this type embeded on its surface has an infinite number 
of members, with their linear dimension increasing as an integer multiple of 
the smallest one. In the case of the tetrahedron, for instance, the first lattice 
has 12 points — those of a truncated tetrahedron — , the next has 48 points, 
and so on following the general rule that the N th generation of the family has 
12iV 2 lattice points. The important point in this construction is that for all 
the members of a given family the curvature is distributed in the same way 
on the lattice. Actually, in the thermodynamic limit we are dealing with a 
kind of discrete version of an orbifold, that is, a degenerate manifold in which 
all the curvature is concentrated at a finite number of points — the vertices 
of the regular polyhedron. When dealing with a given statistical model on 
these lattices, it becomes plausible the physical picture that changing from a 
generation to another inside the same family amounts to consider the same 
system, just with a different length scale. The hypothesis of scaling is suscep- 
tible of being checked on suitable observables, and for this purpose precise 
measurements are presented below in the case of a simple system. 

We focus in what follows on the Ising model formulated over the family 
of lattices embeded on the tetrahedron. In the case of the Ising model we 



2 



may take advantage of powerful techniques developed for the computation of 
partition functions. We refer to the dimer formulation of the problem, which 
we now briefly review || [TD[. Given a collection of spins {<7j}, with i running 



over all the lattice points, the knowledge of the partition function Z requires 
performing the sum over all posible configurations 

Z= e l 3 ^<i,j> a i (X 3 (3) 

o»=±l 

where the sum in the hamiltonian runs over all the lattice links < i,j >. 
The consideration of the high-temperature expansion about (3 = leads to 
the alternative expression 

Z = (cosh J2(tanh [3) n > (4) 
{k} 

where I is the total number of links, stands for the collection of all the 
closed nonintersecting loops on the lattice and {ni} is the respective number 
of links for each of them. The dimer approach manages to compute this sum 
by translating it into the evaluation of Pfaffians of suitable antisymmetric 
operators. This is done as follows. One has first to consider the decorated 
lattice, which in our case is formed by inserting a triangle in place of each 
of the points of the original lattice. A typical decorated lattice looks as that 
depicted in Figure 2, which corresponds to the lattice of the second gener- 
ation. Let us recall the definition of elementary polygon and the notion of 
clockwise odd polygon in the decorated lattice. Elementary polygons are 
those which do not contain lattice points in their interior. An orientation 
may be assigned to each of the links of the decorated lattice, which is repre- 
sented for practical purposes by drawing corresponding arrows on them. A 
polygon is said to be clockwise odd when the number of arrows pointing in 
the clockwise direction on the polygon is odd. Obviously, there is always a 
manner of disposing the arrows on the lattice as to render all the elementary 
polygons clockwise odd. On a planar lattice all the work is done by finding 
one of these systems of arrows. Kasteleyn theorem []11] assures then that the 



combinatorics of the sum in (fD can be reproduced by computing the Pfaffian 
of a single antisymmetric operator A. This is given by a coordination matrix 
on the decorated lattice, such that Aij is positive if an arrow goes from i to 
j and negative in the opposite case. The /3-dependence of the sum is fixed 



3 



by taking the absolute value of matrix elements for triangle links equal to 
z = tank (3 (for /3 > 0), and the absolute value of the rest equal to 1. 

Our lattices embedded on the tetrahedron fall into the category of planar 
lattices, so that Kasteleyn theorem is at work. There are some minor de- 
tails concerning the possibility of defining an order relation and the so-called 
transition cycles on our lattices, which can be elaborated without much com- 
plication. The sum of interest in @ can be represented by the square root of 
the determinant of some antisymmetric operator A following the steps given 
above. There is, in fact, a systematic way of forming a system of arrows 
making all elementary polygons clockwise odd, generation after generation. 
Figure 2 illustrates a particularly convenient way of superposing the deco- 
rated lattices on the plane. The arrows on the triangle links may all be chosen 
clockwise, while those connecting triangles follow the regular pattern shown 
in the bulk. There are only a reduced number of arrows along the boundary 
which cannot conform to this rule and have to be fixed to make all elemen- 
tary polygons clockwise odd. In general, progressing in the family of lattices 
from one generation to the next amounts to add a row of hexagons at the 
top and another at the bottom in the representation of Figure 2, expanding 
accordingly the horizontal dimension. The proposed system of arrows along 
the boundary follows also a regular pattern beginning at each three-fold ring, 
which makes the construction very simple for arbitrarily large lattices. 

Once we have determined the system of arrows for a given lattice, it 
only remains the technical problem of computing the determinant of the 
corresponding antisymmetric operator A. The partition function is given 
now by 

Z = (coshp) 1 (det A) 1 ' 2 (5) 

A direct symbolic computation of the determinant of A cannot be achieved 
in a reasonable amount of time beyond the first generation, which already 
requires to handle a determinant of order 36. We have deviced a trick, how- 
ever, which makes possible the exact computation of the partition function 
at least up to the third generation. The point is that for the kind of operators 
we are dealing with (det A) 1 / 2 is a polynomial in the variable z = tank ft, 
whose determination demands only the knowledge of a finite number of coef- 
ficients. For any of our lattices the degree of the polynomial equals the total 
number of lattice points v. This can be checked by noticing that it always 
exists a closed nonintersecting path which goes through all the points in the 
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lattice. Once it is known that the polynomial (detA) 1 ^ 2 has degree v for each 
lattice, we only need to perform a number of evaluations v of this quantity at 
different values of z to identify all the coefficients (the zeroth order term of 
the polynomial is obviously equal to 1). From a technical point of view, the 
coefficients are integer numbers and can be computed without error by suf- 
ficiently precise evaluations of (det A) 1 / 2 . We have carried out this program 
up to the third generation of lattices by using Mathematica™ . In practice, 
what we have done at each generation is to compute with full precision the 
square root of the determinant at values z — 1, 2, . . . v, for which the result 
of the operation is itself an integer. Then by solving a linear system the 
exact coefficients of the polynomial are obtained. The partition function of 
the Ising model for the first generation lattice of 12 points is, for instance, 

(1 + 4/ + 10/ + 12/ + 15/ + 24/ + 30/° + 24/ 1 + 8/ 2 ) (6) 



1-/ 



For the lattice of the second generation the partition function turns out to 
be 



/ 1 \ 36 

Z = f- — —J (l + 4/ + 28/ + 12/ + 96/ + 96/° + 60/ 1 + 356/ 2 

+576/ 3 + 444/ 4 + 1396/ 5 + 2814/ 6 + 3456/ 7 + 6488/ 8 
+ 12372/ 9 + 18492/° + 29216/ 1 + 50028/ 2 + 77292/ 3 
+115572/ 4 + 177888/ 5 + 265680/ 6 + 383812/ 7 + 545436/ 8 
+760704/ 9 + 1032180/° + 1356444/ 1 + 1728729/ 2 + 2133696/ 3 
+2532312/ 4 + 2870712/ 5 + 3112708/ 6 + 3218976/ 7 + 3142680/ 8 
+2875064/ 9 + 2441460/° + 18981 12/ 1 + 1322784/ 2 + 799344/ 3 
+404928/ 4 + 165888/ 5 + 52272/ 6 + 12096/ 7 + 1728/ 8 ) (7) 

For these first two generations the alternative method of counting polygonal 
paths on the lattice can be worked out, producing the same respective high- 
temperature series as shown above. Our computational approach appears to 
be much more efficient, however, starting with the lattice of the third gener- 
ation. After about 98 hours of CPU time on a Silicon Graphics 4D/480S, we 
have obtained the coefficients of the high-temperature series for the lattice 
of 108 points. The expression for the partition function is given in the Ap- 
pendix. The significance of having at our disposal the analytic expressions 
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of the partition functions is that they allow us to address some questions 
with great precision. One of them refers to the location of the zeroes of 
the partition functions in the z complex plane. These are plotted, for the 
lattices of the second and third generations, in Figures 3 and 4 respectively. 
The frustration of the lattice clearly reflects in the lack of symmetry of both 
partition functions under the exchange (3 —f3. However, the evolution of 
the zeroes in the complex plane points to the formation of two accumula- 
tion points on the real axis in the limit of large lattices, which signal the 
development of two critical points respectively in the ferromagnetic and the 
antiferromagnetic domain. 

The inspection of the first derivatives of the free energy furnishes also clear 
evidence of critical behaviour in the thermodynamic limit. In Figure 5 we 
show a plot of the second derivative of the free energy with respect to f3, for 
the first three generations. At each generation we are able to define a certain 
"critical" coupling constant, within each domain, as that at which the specific 
heat reaches its maximum. Let us focus, for instance, in the ferromagnetic 
domain (/3 > 0), where the critical behaviour has an early development. The 
critical coupling constants for lattices up to the sixth generation (432 points) 
are given in Table 1, along with the corresponding values of the specific heat 
at the maximum. The values for the fourth, fifth and sixth generations result 
from extensive numerical computations carried on the representation ([|) of 
the partition function. The analysis of finite size effects allows to check if the 
scaling hypothesis applies for the lattices in the curved space, being therefore 
helpful in the computation of critical exponents. Regarding the behaviour 
of the finite size critical coupling constant Pl, the standard argument says 
that its difference with respect to the value in the thermodynamic limit (3^ is 
fixed by the point at which the correlation length £ reaches the characteristic 
length dimension L of the lattice. In terms of the v critical exponent we 
should have 

\Pl-Poo\ (8) 

We give in Figure 6 the logarithmic plot of fti — (3^ versus the generation 
number N, that is also the appropriate length scale of the lattice. We have 
assumed that coincides with the critical coupling constant of the Ising 
model on a planar honeycomb lattice, whose value is log(2 + y / 3)/2 ~ 0.6585. 
Quite remarkably, the correlation of the points in the plot is manifest, includ- 
ing those corresponding to the smaller lattices. This provides full support 
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to the scaling law (||). The linear fit for the four bigger lattices in the plot 
gives a value \jv « 1.8, with an uncertainty of about 10%. A significant 
drift from this value should not be discarded for bigger lattices, though the 
investigation of this point requires finer computational techniques. A similar 
analysis can be applied to the divergence of the specific heat. Under the 
assumption of scaling, the maximum value of the finite size specific heat Cl 



has to follow the law [12 



C L ~ L~ a/u (9) 

A logarithmic plot of the values for the specific heat in Table 1 shows again 
a proper linear correlation including the first six generations. The exponent 
which arises from the linear fit is, however, rather small (~ 0.3 taking into 
account the last four generations). The small number of lattices under con- 
sideration makes impossible to discern whether the correct critical behaviour 
follows a power law of the type (|^) or is, in fact, logarithmic. The progression 
to fairly big lattices may be the only way to shed light on this question. 

To summarize, we have proposed an approach to the introduction of cur- 
vature on two-dimensional statistical systems in which the size of the spatial 
lattice may be increased keeping fixed the large scale geometry. We have 
developed it studying the Ising model in particular, showing that the hy- 
pothesis of scaling holds near the critical point. This supports the point 
of view that the thermodynamic limit taken along our sequence of lattices 
provides a sensible definition of the continuum limit of the Ising model on 
a two-dimensional curved space. It would be nice to have an alternative 
field theory description of this limit, specially if one were to confirm unusual 
critical exponents for the specific heat, the spontaneous magnetization or 
the susceptibility of the model. The precise determination of these requires 
further work, as well as a different computational approach. The study of 
finite size effects through the family of lattices proposed appears to be worth 
not only at criticality, but also to investigate possible integrability properties 
off-criticality. 
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tion at the early stages of this work. We want to thank also E. Marino, 
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Appendix 

After the computation described in the text, the partition function for the 
Ising model in the lattice of 108 points turns out to be 



+60Z 11 + lQOlz 12 + 129Qz u + 804z 14 + 5896z 15 + 10U0z lb + 777Qz 17 
+32172^ 18 + 5172(k 19 + 59475^° + 138960^ 21 + 287046^ 22 + 375708^ 23 
+718858^ 24 + 1381752^ 25 + 2154546^ 26 + 3696388^ 27 + 6910026^ 28 
+ 11503728/ 9 + 19206410* 30 + 34305804/ 1 + 58642887^ 32 + 9863796(k 33 
+ 169526346^ 34 + 290348280^ 35 + 486916719^ 36 + 822514464^ 37 + 1387286346^ 38 
+2314063700^ 39 + 3843366849z 40 + 637306956(k 41 + 10492090874z 42 
+ 17155090536^ 43 + 2793596283(k 44 + 45201023376^ 45 + 72645985740/ 6 
+ 116018542776^ 47 + 184128512464^ 48 + 290116289952/ 9 + 453866505468^ 50 
+704934164280^ 51 + 1086527017686z 52 + 1661488438896^ 53 + 2520032972416^ 54 
+3790362808512/ 5 + 565111293324(k 56 + 8349602813680^ 57 + 12221741958408^' 
+ 17717112551232^ 59 + 25427112163896/ + 36115134719952/ 1 
+50748717022488z 62 + 70522684237440z 63 + 96880972147407z 64 
+ 131514489830256/ 5 + 176340456853560/ 6 + 233441767931004/ 7 
+30496303 1957238/ 8 + 392959853131920/ 9 + 499172664304386^ 70 
+624768901981884/ 1 + 770015161230510/ 2 + 933949132164936/ 3 
+ 1114059041938386/ 4 + 1306018538801380/ 5 + 1503582518119761/ 6 
+ 1698618381066144z 77 + 1881426714304464/ 8 + 2041294948806144/ 9 
+2167319364960132/° + 2249464698097728/ 1 + 2279669302012248z 82 
+2252977033783488/ 3 + 2168397154335735/ 4 + 2029388221287936^ 85 
+ 1843829681848794/ 6 + 1623357219696420/ 7 + 1382202439093638^ 88 
+ 1135591499049048/ 9 + 898008996015174/° + 681607478672604/ 1 
+495011279161350/ 2 + 342753592928640/ 3 + 225365558696478/ 4 
+ 140070377148084/ 5 + 81866387456505/ 6 + 44729954075592/ 7 
+22694312231550/ 8 + 10610393197704z 99 + 4531053196539/ 00 
+ 1749121091856/ 01 + 602700352218/ 02 + 182520341604/ 03 
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+47584638729^ 1U4 + 10387733544z 1U5 + 1824896250^ 
+237880152^ 107 + 21484952^ 108 ) 
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N 


Pl 


c L 


1 


0.4673(1) 


1.4392(1) 


2 


0.6082(1) 


1.6520(1) 


3 


0.6342(1) 


1.8924(1) 


4 


0.6439(1) 


2.0798(1) 


5 


0.6486(1) 


2.2308(1) 


6 


0.6513(1) 


2.3568(1) 



Table 1: Respective values of the coupling constant (3l at which the specific 
heat reaches the maximum values Cl for the first six generations. 
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Figure 1: Generic triangular block for honeycomb lattices embeded on the 
tetrahedron, the octaedron and the icosahedron. 
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Figure 2: Decorated lattice for the second generation. The outer lines show 
the identifications of boundary links which embed the lattice on the tetrahe- 
dron. 14 
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Figure 3: Plot in the z-complex plane of the zeros of the partition function 
for the second generation honeycomb lattice on the tetrahedron. 
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Figure 4: Plot in the z-complex plane of the zeros of the partition function 
for the third generation honeycomb lattice on the tetrahedron. 
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Figure 5: Representation of the second derivative of the free energy with 
respect to (5 versus z = tank (3 for the lattices of 12 points (long-dash line), 
48 points (short-dash line) and 108 joints (full line). 




Figure 6: Logarithmic plot of (3l — A» versus the linear dimension L for 
the lattices of the first six generations. The straight line corresponds to the 
linear fit for the last four points. i n 



